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Abstract 

Nineteen  weather  factors  were  analyzed  to  determine  which  have  the 
highest  correlation  with  avalanche  activity.  A  simple  two-parameter  storm  index 
and  a  discriminant  function  model  were  found  to  predict  the  likelihood  of 
avalanches  in  Colorado's  Front  Range.  The  storm  index  utilizes  precipitation 
intensity  modified  by  windspeed  to  predict  the  number  of  avalanches  expected 
on  23  paths.  The  discriminant  function  model  was  used  in  conjunction  with  a 
multiple  regression  program  to  classify  snow  as  stable  or  unstable  on  eight 
avalanche  paths.  This  analysis  shows  that  a  linear  combination  of  the 
maximum  3-hour  precipitation  intensities,  windspeed  resolved  to  an  optimum 
direction  for  each  path,  and  the  sum  of  the  negative  temperature  departures 
from  20°  F  could  be  used  to  determine  whether  a  slope  will  avalanche  70  to  80 
percent  of  the  time. 

Oxford:  384.1:423.5.  Keywords:  Avalanches,  weather,  snow,  statistical 
analysis,  statistical  methods. 


ALL  PHOTOS, 

SEVEN  SISTER  AVALANCHE  PATHS, 
MAY  3,  1973: 

This  slab  avalanche  released  at  4  a.m., 
following  intense  snowfall,  strong  winds, 
and  cold  temperatures.    All  seven  of  the 
north-facing  Seven  Sister  paths  discharged 
their  snow  cover  simultaneously . 

The  fracture  line,  estimated  as  10  feet 
high  in  the  smooth  curved  section,  extended 
across  gullies  and  straight  slopes  for  half 
a  mile  beyond  the  area  shown  at  right  edge 
of  the  centerfold  photo. 

The  tongue-shaped  snow  layer,  barely 
discernible  left  of  the  debris  blocks  (near 
center  of  centerfold  photo),  probably  was 
the  bed  surface  of  a  previous  avalanche. 

The  snow  that  ran  down  the  No.  6 
2 , 000-foot  track  was  12  feet  deep  on  U.S. 
Highway  6. 

Traffic  was  minimal,  and  no  injuries 
or  damage  resulted . 

An  event  of  this  magnitude  rarely 
occurs  in  the  Seven  Sister  area. 

The  use  of  trade  and  company  names  is  for  the 
benefit  of  the  reader;  such  use  does  not  constitute  an 
official  endorsement  or  approval  of  any  service  or  prod- 
uct by  the  U.  S.  Department  of  Agriculture  to  the 
exclusion  of  others  that  may  be  suitable. 
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Predicting  Avalanche  Intensity  from  Weather  Data: 
A  Statistical  Analysis 

Arthur  Judson  and  Bernard  J.  Erickson 


Avalanche  personnel  know  that  certain 
weather  factors  contribute  to  avalanching,  but 
determining  the  relative  contribution  of  various 
factors  is  difficult  for  several  reasons.  For  one 
thing,  most  statistical  analyses  are  complicated 
by  the  high  correlations  between  factors.  Also, 
good  weather  and  avalanche  data  are  rare 
because  they  must  be  continuously  monitored 
in  a  severe  environment,  and  the  complexity 
of  the  relationship  between  weather  and 
avalanches  requires  a  long,  continuous  record. 
Finally,  even  though  avalanche  researchers 
agree  that  precipitation,  wind,  and  temperature 
influence  snow  stability,  they  do  not  agree  on 
the  best  methods  of  quantifying  these  and  other 
factors  for  forecasting  purposes.  The  objectives 
of  this  study  were  to  select  weather  factors  that 
dominate  avalanche  response,  evaluate  their 
relative  contribution,  and  develop  a  predictive 
model. 

Weather  data  for  the  study  were  collected 
at  Berthoud  Pass  in  Colorado's  Front  Range 
from  1951  to  1972.  Avalanche  records  were 
collected  for  the  same  time  period.  Weather 
records  were  taken  by  the  Forest  Service; 
avalanche  data  were  collected  by  the  Forest 
Service,  Colorado  Department  of  Highways,  and 
American  Metals  Urad-Henderson  Mines.  The 
following  weather  data  were  put  on  magnetic 
tape  to  facilitate  analysis  and  to  insure  a  con- 
venient permanent  file:  maximum  3-hour  pre- 
cipitation intensity,  6-hour  average  windspeed 
and  direction,  temperature  and  temperature 
trends,  24-hour  new  snowfall,  water  equivalent, 
total  snow  depth  on  the  ground,  and  tempera- 
ture extremes.  Data  format  for  both  weather 
and  avalanches  is  identical  to  that  used  by 
stations  on  the  Forest  Service  weather  and 
I  avalanche  reporting  network  (Judson  1970). 
'  For  a  first  approximation,  weather  and 
avalanche  data  from  23  avalanche  paths  were 
analyzed  by  univariate  techniques  similar  to 
Perla's  (1970)  analysis  of  contributing  factors 


in  avalanche  hazard  evaluation.  This  univariate 
analysis  resulted  in  a  relatively  simple  two- 
parameter  storm  index  that  effectively  predicts 
the  number  of  avalanches  expected  on  the  23 
paths.  Thirteen  factors  isolated  in  the  univariate 
analysis  were  subsequently  subjected  to  a 
multivariate  discriminant  function  analysis  to 
produce  a  more  refined  three-variable  model  to 
predict  the  likelihood  of  avalanches  on  eight 
paths.  This  model  will  be  thoroughly  evaluated 
and  further  refined  in  the  coming  avalanche 
season. 


The  Univariate  Analysis 


Seven  winters  of  data  (1963-70)  were  used 
in  this  phase  of  the  study.  Twenty-three 
avalanche  paths  located  in  the  Central  Rockies 
near  Berthoud  Pass,  the  Urad  Mine,  and  Love- 
land  Pass  were  selected  for  analysis  (table  1). 
Avalanche  records  for  these  paths  were  uni- 
formly good,  and  the  group  had  terrain  features 
representative  of  most  paths  in  the  Front  Range. 
Nineteen  paths  in  this  group  were  controlled 
by  explosives.  The  group  included  moderate- 
and  high-frequency  paths  which  produced  an 
average  total  of  93  avalanches  per  winter  during 
the  8  winters,  1963-71. 

Scatter  diagrams  and  linear  regression 
proved  useful  in  delineating  important  weather 
factors.  The  number  of  avalanches  from  the 
23  paths  was  plotted  as  a  function  of  single 
weather  factors  or  simple  combinations  of  them 
during  42  storms.  Storms  were  defined  on  the 
basis  of  precipitation  episodes  during  the 
months  of  December  through  March.  The  be- 
ginning of  a  storm  was  defined  as  a  time  when 
measurable  precipitation  fell  in  two  consecutive 
6-hour  periods.  A  storm  was  considered  over 
when  precipitation  ceased  for  three  or  more 
consecutive  6-hour  periods. 
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Table  l.--Some  features  of  the  23  avalanche  paths  (19  controlled)  used  in  the  univariate  analysis, 

7  winters  of  data  (1963-70) 


Path  name 

Locat  ion 

Starting 
zone 
aspect 

Ver t  i  ca 1 
d  rop 

Events  recorded,  1963-71 

Path  Controlled 
frequency  during 
per  winter  period 

Degrees 

Feet 

Number 

Percent 

Lift  Gul ly 

Berthoud  Pass  Ski 

Area 

100 

360 

16 

q4 

Seven  Sister  1 

U.S.  6 

05 

600 

6 

56 

Seven  Sister  3 

U.S.  6 

360 

640 

7 

64 

Seven  Sister  7 

U.S.  6 

360 

750 

2 

81 

Four  B 

Urad  Mine 

105 

1 ,840 

4 

89 

Northwest  Red 

Urad  Mine 

350 

1 ,840 

5 

74 

Rol  1 

Berthoud  Pass  Ski 

Area 

80 

350 

4 

66 

Five  A 

Urad  Mine 

125 

1 ,600 

4 

68 

Five  B 

Urad  Mine 

165 

1 ,810 

3 

75 

Current  Creek 

Berthoud  Pass  Ski 

Area 

85 

500 

8 

0 

Stan  1 ey 

U.S.  ko 

160 

2,760 

3 

54 

Four  A 

Urad  Mine 

90 

1 ,840 

4 

90 

South  Chute  Rol 1 

Berthoud  Pass  Ski 

Area 

30 

350 

3 

74 

One  C 

Urad  Mine 

160 

2,120 

2 

72 

One  E 

Urad  Mine 

170 

1 ,900 

4 

81 

Bethel 

U.S.  6 

160 

2,200 

3 

14 

Dam 

U.S.  40 

70 

2,600 

3 

0 

Five  C 

Urad  Mine 

360 

1 ,840 

2 

71 

Five  Car 

U.S.  6 

75 

400 

2 

33 

Little  Professor 

U.S.  6 

140 

1,360 

2 

27 

Floral  Park 

U.S.  40 

290 

920 

2 

8 

Black  Widow 

U.S.  6 

160 

1 ,640 

2 

17 

Berthoud  Fal 1 s 

U.S.  40 

05 

2,440 

2 

0 

Some  weather  factors,  such  as  total  water 
equivalent  (fig.  1),  showed  definite  trends  while 
others  showed  little  or  no  relation  to  avalanche 
activity  (fig.  2).  Several  factors  were  eliminated 
on  this  basis.  Weather  factors  tested  against 
avalanche  occurrence  by  scatter  diagrams  and 
linear  regression  during  42  storms  were: 

Factors  kept  for  further  analysis 

1.  24-hour  water  equivalent. 

2.  24-hour  snowfall. 

3.  Maximum  precipitation  intensity. 

4.  Maximum  precipitation  intensity  modified 
for  excessive  wind. 

Factors  rejected  from  further  analysis 

1.  Average  windspeed. 

2.  Sums  and  cross  products  of  maximum  pre- 
cipitation intensity  and  windspeed  (several 
combinations). 

3.  Temperature  change  during  storms. 

4.  New  snow  density. 

5.  Settlement. 


Wind  direction  was  not  analyzed  during 
this  part  of  the  study  because  the  paths  had 
many  aspects,  and  wind  direction  varied  con- 
siderably during  single  storms.  Much  of  this 
variation  was  due  to  the  approach  and  passage 
of  upper  level  troughs.  The  usual  sequence  was 
from  SW  to  W  to  NW. 


The  Storm  Index 


Factors  kept  for  further  analysis  were  sub- 
jected to  regression  analysis  using  data  from 
81  storms  during  1963-70.  The  factor  best  cor- 
related with  avalanche  activity  was  the  sum  of 
the  maximum  precipitation  intensities  multi- 
plied by  a  constant  for  excessive  windspeed. 
This  factor  (SPk),  termed  the  storm  index, 
predicts  the  total  number  of  avalanches  (con- 
trolled or  natural)  expected  on  the  23  paths  as 
the  result  of  a  storm.  For  existing  data,  the 
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Figure  1. — Number  of  avalanches  from  23  paths  Figure  2. — Number  of  avalanches  from  23  paths 

as  a  function  of  the  24-hour  water  equivalent  as  a  function  of  the  temperature  change 

of  newly  fallen  snow,  7  winters,  1963-70.  during  storms,  7  winters,  1963-70. 


10  15 
Storm  index 

Figure  3. — The  storm  index,  with  95  percent 
continuous  confidence  intervals  for  the 
mean  of  all  future  observations .  The 
number  of  avalanches  expected,  based  on  7 
winters'  data,  1963-70,  on  the  23  paths  = 

Y  =  -1.31  +  0.76  [e"  Pk] 

where : 

n  =  number  of  6-hour  periods  in  a  storm; 

P  =  maximum  3-hour  precipitation  intensity 

within  each  6-hour  period; 
k  =  a  constant. 

k  =  1  with  windspeeds  <  27  m.p.h. 
k  =  0.3  with  windspeeds  >_  27  m.p.h. 

The  correlation  coefficient  r  =  0.86. 


index  has  a  correlation  r  of  0.86  and  a  standard 
error  of  ±  2.7  (fig.  3). 

The  storm  index  is  computed  every  6  hours 
and  requires  minimum  instrumentation  — a 
recording  precipitation  gage  and  a  continuous 
windspeed  record  give  the  required  data.  This 
index  was  developed  with  avalanche  data  from 
paths  that  have  been  reliable  indicators  of 
avalanche  activity  in  the  Loveland-Berthoud 
Pass  areas  west  of  Denver. 

The  storm  index  was  tested  on  20  storms 
during  1971-72  (fig.  4)  with  satisfactory  results. 


10 

Predicted 

Figure  4. — Scatter  diagram  of  observed  versus 
predicted  avalanches  for  20  storms  during 
1971-72.     Predictions  based  on  the  storm 
index  derived  from  1963-70  data. 
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The  predicted  numbers  of  avalanches  fell  within 
the  confidence  intervals  shown  in  figure  3. 
Observed  values  were  not  consistently  above 
or  below  the  predicted,  and  scatter  about  the 
regression  line  appears  to  be  random.  Much  of 
the  scatter  is  attributed  to  nonuniform  control 
efforts  due  to  other  operational  considerations 
confronted  by  highway  and  mine  personnel. 

It  was  interesting  to  note  that  another 
storm  index,  developed  with  avalanche  data 
from  23  uncontrolled  avalanche  paths,  had  a 
much  lower  correlation  coefficient  and  a  greater 
degree  of  scatter  than  the  one  we  finally 
developed.  This  difference  implies  that  fore- 
casting indices  based  on  avalanche  data  from 
uncontrolled  paths  are  difficult  to  interpret  and 
are  less  reliable  as  forecast  guides. 

The  main  drawback  with  the  storm  index 
is  that  the  index  is  highest  near  the  end  of 
storms,  even  though  hazard  may  be  decreasing 
because  some  avalanches  have  already  fallen 
and  the  snow  is  stabilizing.  A  way  of  reducing 
the  index  toward  the  end  of  the  storm  (a  decay 
function)  is  badly  needed  and  is  now  under 
study. 


Applications 


The  storm  index  should  have  some  applica- 
tion in  other  mountain  areas.  It  utilizes 
precipitation  intensity  and  windspeed  during 
storms  to  predict  the  expected  intensity  of 
avalanching  on  an  area  basis.  Because  this 
index  was  developed  using  avalanche  data  from 
23  paths  that  have  a  wide  variety  of  physical 
features,  we  believe  it  will  work,  with  some 
modifications,  at  most  avalanche  areas.  The 
two  weather  factors  comprising  the  storm  index 
are  directly  related  to  the  rate  of  loading  on 
avalanche  slopes.  The  rate  of  loading  is  a  prime 
factor  contributing  to  avalanche  release  at  any 
area.  Plans  are  being  developed  for  testing  the 
storm  index  at  other  avalanche  areas. 


The  Multivariate  Analysis 


Weather  factors  other  than  precipitation 
intensity  and  windspeed  affect  avalanche  forma- 
tion. Moreover,  because  it  is  the  combined  effect 
of  several  factors  which  determines  snow 
stability,  it  appeared  logical  to  try  a  multivariate 
approach  to  predicting  avalanche  potential. 
Also,  factors  affect  individual  paths  in  different 


ways,  so  a  single  path  analysis  was  indicated. 
A  discriminant  function  analysis  was  selected 
for  the  second  phase  of  the  study. 

A  group  of  10  well-defined,  controlled  paths, 
(table  2)  were  selected  for  this  analysis.  Data 
from  1952-71  were  analyzed.  Like  those  in  the 
first  phase  of  this  study,  these  paths  run  fre- 
quently, and  data  on  their  occurrence  and 
control  were  uniformly  good.  They  are  repre- 
sentative of  many  moderate-  to  high-frequency 
paths  in  the  Front  Range  because  of  their 
location,  physical  features,  and  wide  variety  of 
starting  zone  aspects.  Eight  paths  threaten 
highways  or  roads,  and  two  are  in  a  ski  area. 
None  are  skied. 


The  Discriminant  Function 


The  discriminant  function  is  a  multivariate 
statistical  technique  of  assigning  data  into  two 
or  more  groups  based  on  prior  knowledge.  More 
specifically,  it  is  the  linear  component  of  p 
variables  which  maximizes  the  ratio  of  the 
between-groups  variance  to  the  within-groups 
variance. 

In  general  form,  the  discriminant  function 
is  written  as 

L  =  PiXi  +  SzXz  +  .    .    .  +  6  X  [1] 

r  It 

The  discriminant  coefficients  3i,  Ba  •  •  •  Sp  are 
computed  using  the  Fisher  method  discussed 
by  Rao  (1952).  The  multivariate  mean  of  group 
1  is  given  by 

Ri  =  BiXii  +  B2X21  +  .   .   .  +  B  f  [2] 

P  Pi 

and  of  group  2  by 

R2  =  B1X12  +  B2X22  +  .    .    .   +  BpXp_^  [3] 

The  generalized  distance  between  the  group 
means,  called  Mahalanobis'  D^,  is  the  difference 
between  the  group  means.  Hence 

=  Ri  -  R2  [4] 

The  discriminant  function  is  tested  for  signif- 
icance using  the  F  ratio  involving  D^.  When 
the  function  is  significant,  there  is  a  real 
difference  between  groups.  The  discriminant 
index  Rq,  which  determines  the  group  classifica- 
tion of  future  data,  is  a  weighted  average  of 
the  group  means: 

niRi  +  n2R2  [5] 
ni  +  n2 
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Table  2. — Pertinent  features  of  the  10  avalanche  paths  (controlled)  used  in  the  multivariate 

analysis,   1952-71  data 


,,     ^.     Events  recorded,  1963-71 

^.  Verti-   ' 

r,  1       ^-  btartinq  zone  ,   ;^ — — r  ? — t — n — T 

Path  name  Location   cal  Path  Controlled 

Area    Aspect  Shape  drop    frequency  during 

per  winter  period 

Acres  Degrees  Feet        Number  Percent 

Lift  Gully  Berthoud  Pass      0.2     100        Bowl-shaped  360         16  Sk 

Ski  Area  depression 

Cliff  Berthoud  Pass      0.^*      90       Straight  ramp  350  9  89 

Ski  Area 

Floral  Park         U.S.  40  10        290        Poorly  defined,  920  2  8 

si ight  depressions 
with  a  midway  bench 

Stanley  U.S.  40  20       160       Broad  bowl-shaped  2,760  3  54 

depressions  and 
shal low  gul 1 ies 

Northwest  Red      Urad  Mine  22       350       Steep  bowl-shaped  1  ,840  5  74 

depression  with 
gul 1 ies 

Four  A  Urad  Mine  10         90       Shallow  gul 1 ies  1,840  4  90 

topped  by  a  cl i f f 

Four  B  Urad  Mine  15        105       Bowl-shaped  1 ,840  4  89 

depression  with 
centra  1  gully 

Bethel  U.S.  6  15        160       Straight  slope  2,200  3  14 

feeding  a  gul ly 
from  the  side 

Seven  Sister  3    U.S.  6  0.5    360       Broad  and  shallow  640  8  64 

gul  ly 

Seven  Sister  6    U.S.  6  4       360       Sha 1  low  depress  ion         1,050  7  56 

flanked  by  a 
prominent  rock  rib 


When  L  >  Rq,  the  event  is  classified  in  group  1. 
The  vahdity  of  Rq  is  attained  by  computing  the 
probability  of  misclassification.  This  is  done  by 
entering  D/2  in  a  cumulative  normal  frequency 
distribution  table  of  the  normal  deviate. 

The  discriminant  function  was  introduced 
by  Fisher  in  1938.  It  has  since  been  used  in 
medicine,  psychology,  engineering,  biology, 
economics,  anthropology,  and  geology.  Weeks^ 
was  first  to  use  the  technique  on  weather  and 
avalanche  data.  Bois  and  Obled  [1972]  recently 
began  work  with  the  discriminant  function  on 
data  from  Switzerland. 


''Unpublished  data  [1967],  U.  S.  Army  Materiel 
Command,  Cold  Regions  Res.  and  Eng.  Lab.,  Hanover, 
N.H. 


Analysis  of  Data 

This  study  was  confined  to  the  period, 
November  15  to  April  15,  which  limits  the 
analysis  to  dry  snow  at  this  site.  A  20-inch 
snow  depth  threshold  at  the  Berthoud  Pass 
weather  and  snow  study  site  was  required 
before  analysis  began.  Data  from  1952  through 
1971  were  used.  Group  classifications  were  based 
on  control  results.  Snow  was  defined  as  un- 
stable (group  1)  when  control  efforts  produced 
a  slide  or  when  a  natural  avalanche  occurred. 
Snow  was  classified  stable  (group  2)  when 
control  efforts  failed  to  initiate  an  avalanche. 
Weather  data  for  both  groups  were  limited  to 
the  interval  between  control  applications  and/or 
to  the  time  between  natural  avalanches  on  each 
path.  These  intervals  varied  from  1  day  to  6 
weeks. 


5 


A  preprocessing  program  (WXDAT)  was 
written  to  select  weather  data  for  both  groups 
from  a  master  file  of  Berthoud  Pass  raw 
weather  data  on  magnetic  tape.  WXDAT  pre- 
pares an  input  file  of  selected  weather  variables 
for  the  discriminant  function  program.  The 
program  written  by  Davis  and  Sampson  (1966) 
for  the  IBM  1620  was  used  in  our  analysis.  This 
program  was  converted  for  use  on  a  CDC  6400 
computer  at  Colorado  State  University.  The 
program  avoids  complex  matrix  inversion  and 
computes  the  discriminant  function  coefficients 
for  two  groups  with  a  maximum  of  20  variables. 
Program  outputs  include  the  discriminant  index 
Rq,  the  multivariate  group  means  R,  and  R2, 
the  F  ratio,  and  Mahalanobis'  D^.  Additional 
outputs  were  added  as  analytical  aids.  The 
number  of  cases  in  each  group  does  not  have 
to  be  equal,  but  must  exceed  the  number  of 
variables. 

Variates  used  with  the  program  must  not 
be  highly  interrelated,  a  constraint  which 
presented  an  immediate  problem  since  weather 
variables  are  correlated.  Fortunately,  a  linear 
stepwise  program  could  be  applied  to  the  dis- 
criminant problem  due  to  the  mathematical 
equivalence  between  the  models.  The  correla- 
tion matrix  provided  by  the  stepwise  program 
was  used  to  delete  highly  interrelated  variates 
which  add  little  to  the  analysis. 

For  computational  convenience,  the  step- 
wise program  was  given  a  dummied  predictand 
as  follows:  Taking  n,  cases  in  group  1  and  n2 
cases  in  group  2,  we  compute  a  predictand 

/  n,  4-  n2  for  all  cases  in  group  1  and 
-n,  /  n,  -H  n2  for  all  cases  in  group  2.  The 
average  value  of  the  predictand  is  zero  since 


n2 

-1-  nj 

-n. 

n,  -1-  n2 

n,  -1-  n2 

=  0 


[6] 


Used  in  this  manner,  the  stepwise  or  screening 
program  provides  a  ranking  of  variates  based 
on  the  increase  in  explained  variance.  The  first 
variate  is  selected  on  the  basis  of  having  the 
highest  F  value  at  a  given  level.  Once  selected, 
it  is  held  aside,  and  the  program  selects  the 
variate  with  the  highest  F  value  of  the 
remaining  parameters,  and  so  on.  The  screening 
program  provides  an  unbiased  selection  of 
variates  which  are  then  run  through  the  dis- 
criminant function  program.  The  screening 
program  is  terminated  when  the  additional 
variables  entered  are  not  significant. 

The  following  13  variables  were  selected 
for  analysis  on  the  basis  of  field  experience 
and  the  earlier  univariate  analysis: 


1.  Sum  of  24-hour  water  equivalents. 

2.  Sum  of  the  maximum  3-hour  precipitation 
intensities  in  each  6-hour  period  decayed 
over  the  interval. 

3.  Same  as  No.  2  without  the  decay  function. 

4.  Sum  of  maximum  precipitation  intensities 
decreased  for  excessive  windspeeds. 

5.  Count  of  6-hour  windspeeds  <  15  m.p.h. 

6.  Count  of  6-hour  windspeeds  between  15  and 
27  m.p.h. 

7.  Count  of  6-hour  windspeeds  >  27  m.p.h. 

8.  Count  of  6-hour  temperatures  >    20°  F. 

9.  Sum  of  the  negative  temperature  departures 
from  20°  F. 

10.  Average  precipitation  intensity. 

11.  Sum  of  windspeeds  during  precipitation 
periods,  resolved  to  an  optimum  direction 
for  each  path. 

12.  Sum  of  windspeeds  resolved  to  an  optimum 
direction  for  each  path. 

13.  Sum  of  squares  of  windspeeds  resolved  to 
an  optimum  direction  for  each  path. 


Results  and  Discussion 


The  discriminant  function  on  8  of  the  10 
paths  was  significant  at  either  the  1-  or  5-per- 
cent level.  The  avalanche  records  for  the  two 
paths  failing  the  significance  test  were  found 
to  be  incomplete  because  control  teams  failed 
to  enter  negative  control  results.  Such  data  are 
critical  since  they  define  the  time  interval  for 
weather  data  used  in  the  analysis.  Neither  path 
provided  a  meaningful  discriminant  function. 

Results  of  the  study  are  summarized  in 
table  3.  The  most  important  variates  are: 


No.  Designation 


2PI*D 


9      ZNeg  TT 


12 


2VVr 


Variable 

Sum  of  the  maximum  consec- 
utive 3-hour  precipitation 
intensities  within  each  6-hour 
period  decayed  over  the 
interval. 

Sum  of  the  6-hour  negative 
temperature  departures  from 
20°  F. 

Sum  of  the  windspeeds  >  15 
m.p.h.  resolved  to  an  optimum 
direction  for  each  path. 


Maximum  precipitation  rates  with  a  decay  func- 
tion (2PI*D)  dominated  avalanche  response  on 
five  of  the  eight  paths,  while  windspeed 
resolved  to  an  optimum  direction  for  each  path 
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Table  3- 

--Compos  i 

te  summary 

of  statistical 

data,  mu 

tivar iate  analysis, 

1952-71  data 

Path 
(1 ) 

Signifi- 
cance 
1  evel 
(2) 

Coef f  i  c  i ■ 
6  n  t 

(3) 

Variable 

3  yiM tj\j  I 

(4) 

Var  iabl e 

n  1  im       r"  ^ 
1 1  Ulll  UC  1 

(5) 

(6) 

Probabi 1 i  ty 
of  misclassi- 
f  i cat  ion 
(7) 

Increase  in 
expla  ined 
variance^ 
(8) 

R 

o 

(9)  1 

Percent 

Lift  Gul ly 

1 

0.01 1 

.280 
.000 

ZPhvD 
ZNeg  TT 

12 

2 
9 

1 .68 

0.26 

0.1907 

.0967 
.0001 

2.26  ' 

CI  iff 

1 

.416 
.004 

EPI^D 

2 
12 

2.06 

.24 

.2904 
.0446 

3.37 

Floral  Park 

5 

.495 
.001 
.000 

EP  1  '^D 
ENeg  TT 

2 
12 
9 

2.50 

.21 

.3909 
.0126 
.0021 

4.64 

Stanley 

1 

.294 
.004 
-  .004 

EPh'^D 
ZNeg  TT 

2 
9 
12 

1 .60 

.26 

.1863 
.0640 
.0438 

2.66 

Northwest  Red 

1 

.005 

.141 
.000 

IP  1*0 
ZNeg  TT 

12 

2 
9 

1.07 

.30 

.1671 

.0440 
.0053 

1.61 

Four  A 

5 

.288 
.003 
.000 

ZPI*D 

INeg.TT 

EVV,,P 
w 

2 
9 

12 

1 .06 

.29 

.1655 
.0530 
.0005 

2.73 

Four  B 

1 

.011 

.004 
.084 

ZNeg  TT 
ZPI^D 

12 

9 
2 

2.10 

.24 

.3095 

.0295 
.0119 

2.97 

Bethel 

1 

.325 
.003 

.000 

ZPI^D 

NW 
ZNeg  TT 

2 

12 
9 

2.23 

.23 

.3048 
.0407 

.0141 

4.13 

^Variable  symbols  and  numbers  are  explained  on  page  9. 

^Numbers  in  this  column  were  provided  by  the  screening  program;  they  determine  the  relative 
order  of  importance  of  the  variables. 


(2VVj^)  was  the  primary  factor  on  three  paths. 
An  arbitrary  decay  function  was  used  to  decrease 
the  sum  of  the  precipitation  intensities  with 
time.  The  function  is  held  at  one  for  the  first 
2  days,  reaches  0.5  on  the  5th  day,  and  levels 
off  at  0.2  from  the  9th  day  on.  This  function, 
used  to  simulate  stabilization  with  time, 
decreased  the  probability  of  misclassificationon 
seven  of  eight  paths.  This  parameter  is  currently 
being  refined  for  both  wind  and  precipitation. 
Temperature  appears  to  be  an  important  second- 
ary factor.  It  was  the  second  most  important 
variate  on  three  paths. 

The  stepwise  program  is  affected  by  cor- 
relation between  independent  variables.  There- 


fore, the  order  of  importance  (table  3)  is 
relative.  The  primary  validity  of  the  variables 
selected  is  given  by  the  probability  of  misclass- 
ification  in  column  7  of  table  3.  Variables 
contributing  the  least  toward  the  explained 
variance  have  coefficients  near  zero.  Because 
there  is  some  correlation  between  variables, 
coefficients  may  be  negative  in  a  few  cases, 
even  though  the  physical  effect  of  the  variable 
is  assumed  to  be  positive. 

The  multivariate  analysis  indicates  that  (1) 
precipitation  intensity,  (2)  windspeed  resolved 
to  an  optimum  direction  for  each  path,  and  (3) 
temperature  can  be  used  to  predict  the  likelihood 
of  avalanches. X  Even  more  important,  these  3 
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,  variables  predict  avalanche  activity  on  the  test 
;  paths  more  accurately  than  do  all  13  variables 
:  combined. 

The  variables  and  order  of  importance  are 
I  different  for  different  paths.  Paths  with  a  high 

frequency  of  occurrence  are  less  affected  by  low 

temperatures  than  are  paths  that  run  less  often, 
I  which  is  no  surprise  since  it  takes  time  to 

change  snow  structure.  The  effects  of  rapid 
1  temperature  changes  on  avalanche  release  were 
mot  examined  due  to  insufficient  data.  In  this 

regard,  one  must  realize  that  rapidly  falling 
I  temperatures  occur  almost  every  afternoon, 
I  but  avalanches  do  not. 


Applications 


The  technique  and  analysis  developed  in 
the  second  phase  of  this  study  can  be  applied 
to  weather  and  avalanche  data  anywhere. 
Weather  factors  dominating  avalanche  activity 
at  the  Colorado  study  area  will  probably  be  key 
factors  in  other  areas  where  dry  snow  avalanches 
are  the  main  problem.  The  specific  coefficients 
derived  in  this  study  will  be  different  at  other 
locations,  but  the  magnitude  of  these  differences 
is  not  presently  known. 


Current  plans  call  for  testing  the  three- 
variable  model  by  using  the  eight  avalanche 
paths  used  in  the  second  phase  of  this  study  as 
an  index  of  avalanche  activity  in  the  100-square- 
milearea  between  Arapaho  Basin  and  Berthoud 
Pass.  The  response  from  the  index  group 
resembles  the  occurrence  pattern  of  avalanches 
from  40  other  paths  in  the  area  (fig.  5).  The 
probability  of  misclassification  (column  7,  table 
3)  indicates  activity  from  the  index  group  can 
be  correctly  predicted  70  to  80  percent  of  the 
time.  If  an  independent  set  of  data  confirm 
this  accuracy  the  discriminant  function  coeffi- 
cients and  variables  for  the  index  group  could 
serve  as  a  basis  for  regional  avalanche  warnings 
as  well  as  a  guide  for  avalanche  control  on  the 
eight  index  paths. 

The  L  values,  or  discriminant  scores,  could 
be  calculated  and  pooled  for  the  index  paths 
every  6  hours  to  simulate  avalanche  activity 
on  a  real-time  basis.  When  L  exceeds  the  dis- 
criminant index  R©  for  a  given  path,  the 
function  predicts  the  path  will  avalanche  either 
naturally  or  artificially.  For  example,  the  dis- 
criminant function  for  the  Lift  Gully  in  table  3 
is: 


L  =  O.OlKZWj^^)  +  0.28(ZPI*D). 
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Figure  5. — Relationship  between  daily  avalanche  response  from  the  8  index 
paths  and  40  other  paths  in  the  same  general  area,  March  10  to  April  10, 
1965. 
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During  one  42-hour  period  in  December  1964, 
the  values  for  the  variables  were  65.9  and  8.5, 
respectively.  Substituting  in  the  formula: 

L  =  0.011(65.9)  4-  0.28(8.5)    =  3.1 

Ro  for  this  path  is  2.26,  the  discriminant  index 
was  exceeded,  and  the  slope  avalanched.  Similar 
calculations  can  be  made  for  all  index  paths. 

At  present,  individual  L  values  are  reset  to 
zero  when  a  slope  avalanches  or  when  it  is  shot 
with  no  results.  The  decay  function  reduces 
the  sum^  of  the  maximum  precipitation  intensi- 
ties (2PI)  with  time.  Wind  and  temperature 
could  be  treated  in  the  same  manner,  so  that 
L  for  all  paths  and  therefore  L  -  Rq,  which 
represents  the  likelihood  of  avalanching,  would 
fluctuate  through  periodic  cycles  simulating 
avalanche  activity  in  the  region. 


Summary  and  Implications 


The  storm  index  developed  and  tested  in 
the  first  phase  of  the  study  consists  of  the  sum 
of  the  maximum  3-hour  precipitation  intensities 
multiplied  by  a  constant  for  excessive  wind- 
speed  (2Pk).  It  predicts  the  number  of  ava- 
lanches expected  on  23  paths  in  Colorado's 
Front  Range  during  storms.  It  could  be  used  as 
an  objective  guide  for  issuing  regional  avalanche 
warnings.  Its  main  limitations  are:  (1)  it  can 
be  used  only  during  precipitation  periods,  and 
(2)  it  contains  no  provision  for  a  decrease  in 
avalanche  hazard  during  and  following  storms. 
With  modifications,  it  can  probably  be  used  at 
other  mountain  locations  where  dry  snow  ava- 
lanches are  the  primary  problem.  An  identical 
model  made  with  the  same  weather  data  but 
based  on  avalanche  activity  from  23  uncon- 
trolled paths  yielded  poor  results.  It  is  therefore 
recommended  that  similar  models  developed 
at  other  mountain  areas  be  based  on  avalanche 
activity  from  paths  which  are  controlled. 

The  three-variable  model  developed  in  the 
second  phase  uses  maximum  precipitation  in- 
tensity, windspeed  and  direction,  and  tempera- 
ture to  predict  the  likelihood  of  avalanches  on 
eight  controlled  paths  intheBerthoud-Loveland 
Pass  area.  The  model  is  untested.  If  planned 
test  results  are  satisfactory,  this  model  will 
serve  as  a  guide  for  control  decisions  on  the 
eight  paths.  It  could  also  provide  the  means 
for  issuing  a  regional  avalanche  warning  and 
for  lifting  that  warning.  The  three-variable 
model  is  more  flexible  than  the  storm  index 


because  the  calculation  of  L,  or  the  likelihood 
of  avalanching,  begins  when  snow  depth  on 
the  ground  reaches  a  threshold  depth  in  early 
winter,  and  continues  uninterrupted  until 
spring.  The  L  values  increase  and  decrease  with 
time  and  avalanche  activity.  This  model  could 
provide  a  daily  rating  of  avalanche  hazard 
which  is  independent  of  any  subjective  defini- 
tion of  storms.  The  main  contribution  from  this 
second  phase  is  the  technique  which  can  be 
used  to  evaluate  weather  and  avalanche  data 
from  any  mountain  area.  The  weather  factors 
which  were  found  to  dominate  avalanche  activity 
in  this  study  will  probably  be  key  weather 
factors  affecting  avalanche  response  at  other 
mountain  locations  where  dry  snow  is  common. 

Although  much  work  remains,  an  objective 
basis  for  issuing  regional  avalanche  warnings 
is  now  possible.  Additional  variables  that  would 
give  an  indication  of  stress  conditions  within 
the  snow  cover,  if  they  can  be  isolated  and 
measured,  would  enhance  the  reliability  of  such 
warnings.  Until  we  learn  more  about  how  to 
isolate  and  measure  the  snow-cover  features 
pertinent  to  avalanche  release,  any  objective 
avalanche  warning  scheme  will  be  primarily 
dependent  on  weather  factors. 


Literature  Cited 


Bois,  Philippe,  and  Charles  Obled. 
[1972].  Analyse  des  donnees  nivoclimato- 
logiques  en  vue  de  la  prevision  des  ava- 
lanches. 71  p.  Institut  federal  pour  I'etude 
de  la  neige  et  des  avalanches.  Weissfluhjoch 
-  Davos,  Switzerland. 

Davis,  John  C,  and  Robert  J.  Sampson. 
1966.   Fortran  II  program  for  multivariate 
discriminant  analysis  using  an  IBM  1620 
computer.  Kans.  Geol.  Surv.  Comput.  Con- 
trol 4,  8  p. 

Judson,  A. 

1970.  A  pilot  study  of  weather,  snow,  and 
avalanche  reporting  for  western  United 
States,  p.  123-134.  In  Ice  engineering  and 
avalanche  forecasting  and  control.  Nat.  Res. 
Coun.  Can.  Tech.  Memo.  98,  163  p. 
Perla,  R.  I. 

1970.    On  contributory  factors  in  avalanche 
hazard  evaluation.  Can.  Geotech.  J.  7(4): 
414-419. 
Rao,  C.  R. 

1952.  Advanced  statistical  methods  in  bio- 
metric  research.  390  p.  John  Wiley  and 
Sons,  N.Y. ;  Chapman  and  Hall,  Ltd.,  London. 


Agricukyre— CSU,  Ft.  Collins 


12 


